Structure factor and thermodynamics of rigid dendrimers in solution 
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The "polymer reference interaction site model" (PRISM) integral equation theory is used to 
determine the structure factor of rigid dendrimers in solution. The theory is quite successful in 
reproducing experimental structure factors for various dendrimer concentrations. In addition, the 
structure factor at vanishing scattering vector is calculated via the compressibility equation using 
scaled particle theory and fundamental measure theory. The results as predicted by both theories 
are systematically smaller than the experimental and PRISM data for platelike dendrimers. 
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I. INTRODUCTION 



be embodied in the structure factor S{q, p) defined as 



Dendrimers are synthetic macromolecules with defined 
architecture that are synthesized by iterative controlled 
reaction steps. Starting from a focal unit, subsequent 
generations are connected to this initial core which re- 
sults in a tree-like structureii^i^i^. Today it is well- 
understood that dendrimers composed of flexible units 
adopt a so-called dense core structure, that is, the aver- 
age segmental density has its maximum at the center of 
the molecule^i^i^. This is easily derived from the fact that 
flexible dendrimers can assume a great number of con- 
formations in which the terminal groups can fold back. 
Hence these flexible dendrimers do not exhibit a well- 
deflned surface given by the terminal groups of the last 
generation. The average density proflle thus derived can 
be used to calculate the interaction of flexible dendrimers 
in solution. Thus these structure are well-understood by 
now^ . 

Much less attention has been paid to dendrimers con- 
sisting of rigid units^i^'^iiSiiii'i^iii^'i^. Figured] (a) shows 
the chemical structure of such a rigid dendrimer which is 
solely composed of stiff units^i^iSiiiS. While the dendritic 
scaffold of flexible dendrimers in solution can adopt a 
large number of conformations which follow from rota- 
tions about various bonds, dendrimers consisting of such 
rigid units exhibit a rather well deflned structure in so- 
lution. This fact has been shown recently by small-angle 
neutron scattering (SANS)^'^'^^. Hence, these systems 
may serve as model systems for interacting monodisperse 
particles in statistical physics. 

Here we present a comprehensive discussion of the in- 
teraction of rigid dendrimers in solution. Two systems 
involving rigid dendrimers are examined in the present 
paper, namely polyphenylene dendrimers of the fourth 
generation in solution (see Fig.[T](a)) and stilbenoid den- 
drimers of the third generation in solution (see Fig.[T](b)). 
The problem to be addressed is that of structural proper- 
ties over a range of dendrimer concentrations. In general, 
the effect of mutual interaction of dissolved species with 
number density p = N/V in a scattering experiment can 



^Ioiq)S{q,p) = yV^'iAp)'P{q)Siq,p) (1) 



V 



where Iq (q) is the scattering intensity of the single par- 
ticle, Vp the volume of the particle, P{q) the form factor 
(normalized to unity for q — 0) and q is the magnitude of 
the scattering vector {q — {Att / \)sin{9 /2), A: wavelength 
of radiation, 9: scattering angle). Ap is the contrast of 
the solute resulting from the difference of the average 
scattering length density and the scattering length den- 
sity of the solvent. Up to now, mutual interaction has 
only been addressed because scattering experiments re- 
quire flnite concentrations in order to obtain data with 
a reasonable statisticaii^iii^. Hence, the effect of mutual 
interaction of the dissolved dendrimers has only been con- 
sidered when extrapolating the data to vanishing concen- 
trations. Details of such a data evaluation may be found 
in Refi^. Here we determine the correlation functions 
and structure factors S{q, p) within the framework of an 
interaction site model and compare these results to ex- 
perimental d&ta^^^. Moreover, the validity of the pre- 
dictions of both scaled particle theory and fundamental 
measure theory are investigated. A systematic compari- 
son is made between the results of the theories and ex- 
perimental data for two differently shaped dendrimers in 
solution, in order to quantify the influence of the shape of 
the particles on the quality of the predictions of the theo- 
ries. We demonstrate that rigid dendrimers present a new 
class of strictly monodisperse model colloids whose inter- 
action can be conveniently studied in solution. In this 
way we show that these systems allow us to re-consider 
the basic statistical mechanics of rigid objects in solution 
and to compare theory and experiment for the first time 
in a rigorous fashion. 

The paper is organized as follows: In the next section 
we briefly recapitulate the analytical models necessary to 
calculate the scattering intensity for a system of interact- 
ing rigid bodies. Thereafter the theory will be compared 
to experimental results obtained for the rigid dendrimers 
shown in Figs. [T] (a) and (b). In section IV these results 
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FIG. 1: (a) Chemical structure of the polyphenylene den- 
drimer of the fourth generation, (b) Chemical structure of the 
stilbenoid dendrimer of the third generation (R = OCe-ff 13). 



are compared to the fundamentals of the thermodynam- 
ics of rigid bodies. Special attention is paid to a discus- 
sion of exact virial coefficients. A brief conclusion will 
summarize the results at the end. 



II. RISM AND PRISM THEORIES 

The systems under investigation are solutions. How- 
ever, in view of the mesoscopic scale of the particles, the 
solvent will be considered as a structureless continuum. 
Spatial pair correlations of an isotropic fluid of identi- 
cal particles, each carrying n distinct interaction sites, 
are characterized by a set of intermolecular site-site total 
correlation functions hij{r,p), where the indices i and j 
run over sites on each of two particles and p is the par- 
ticle number density. These functions are related to a 



set of intermolecular site-site direct correlation functions 
Cij (r, p) by the generalized Ornstein-Zernike relations of 
the "reference interaction site model" (RISM), which in 
Fourier space rea d^^'^^ 



(g, p)cmo{q, p) {(^ojiq, p) + phojiq, p)) 



where the (jJij{q,p) are the Fourier transforms of the in- 
tramolecular correlation functions. The set of generalized 
Ornstein-Zernike equations must be supplemented by a 
set of closure relations. If the interaction sites are simply 
the centers of exclusion spheres of diameter d, to account 
for steric effects, a convenient closure is the Percus-Yevick 
approximatio] 



.16.18 



hj{r,p) = -I, r <d, Cij{r,p) = 0, r > d. (2) 

The experimentally accessible structure factor S{q, p) is 
defined as 



Siq,p) = 1 + p 



h{q,p) 



where 



Piq,p)' 
1 " 



(3) 



(4) 



is the particle-averaged total correlation function. The 
particle-averaged intramolecular correlation function 



1 " 

nq.P)--, E 



M^p) 



(5) 



characterizes the geometry of the distribution of the sites, 
and hence the geometric shape of the particles. While 
the particle-averaged intramolecular correlation function 
accounts for the interference of radiation scattered from 
different parts of the same particle in a scattering ex- 
periment, the local order in the fluid is characterized by 
h{q,p) or S{q,p). 

The RISM has been proved to be a successful theory of 
the pair structure of many molecular fluids (for a review 
see Refji^). In the case of macromolecular and colloidal 
systems, with very large numbers of interaction sites, the 
number of coupled RISM equations becomes intractable, 
and a considerable simplification follows from the as- 
sumption that the direct correlation functions Cy {q, p) 
are independent of the indices i and j. This leads to 
the "polymer reference interaction site model" (PRISM) 
theory first applied by Schweizer and Curro to long flexi- 
ble polymers^*^ . PRISM neglects end effects in that case. 
The resulting single generalized Ornstein-Zernike equa- 
tion of the PRISM reads 

Hq, P) = c{q, p)P'^{q, p) + pc{q, p)h{q, p)P{q, p) , (6) 

where c{q, p) = YZi,o=i'^rao{q, p)- The PRISM integral 
equation theory has been successfully applied to vari- 
ous systems, such as rodlike viruses2ii^i^i2i, platelike 
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FIG. 2: The form factor P{q) — P{q, p 0) of polyphenylene 
dendrimers of the fourth generation^ (squares and upper hne) 
and stilbenoid dendrimers of the third generations^ (circles 
and lower line) as obtained by small angle neutron scattering 
(symbols). The lines represent the calculated form factors. 
For reasons of clarity, the upper data set and line have been 
shifted up. 



folds set up from stilbenoid unitgiiii^iii^. Starting from a 
central phenyl group all subsequent generations are built 
up by tra77.s-stilben units, only the terminal groups are 
substituted by hexyloxy groups in order to ensure better 
solubility in common solvents such as toluene (see Fig. [1] 
(b)). Full conjugation in trans-stilbene can be achieved 
in the completely planar conformation. However, the po- 
tential energy for a slight torsion around the single bonds 
is low in the ground state of trans-stilbene. Both molecu- 
lar modelingi^ and a small-angle scattering studyi^ have 
demonstrated that a stilbenoid dendrimer of the third 
generation exhibits a relatively compact platelike struc- 
ture. 

The measured and calculated form factors P{q) of both 
polyphenylene dendrimers of the fourth generation and 
stilbenoid dendrimers of the third generation are shown 
in Fig. [21 The form factor of the stilbenoid dendrimers 
agrees with the one calculated numerically for a cir- 
cular platelet with a radius R = 2.4 nm and thickness 
L — 1.& nm according to 



colloids^^i^i^ and dendrimers^'*, polymers22i2^, mixtures 
of spherical colloids and semiflexible polymers^, and 
bottlebrush polymers'^"'^ . Moreover, it has been demon- 
strated recently that the simpler PRISM theory yields 
results in good agreement with the more elaborate RISM 
calculations for lamellar coUoids^^. 

For flexible polymers or dendrimers the particle- 
averaged intramolecular correlation function P(q, p) de- 
pends on the particle number density and follows from 
a statistcal average over particle configurations. In 
the limit p — > the particle-averaged intramolec- 
ular correlation function reduces to the form factor 
P{q) = P{q, p 0). In the case of rigid particles P{q, p) 
is independent of the particle number density because the 
particles are not deformed due to intermolecular interac- 
tions for typical concentrations in the fluid state. 



III. RIGID POLYPHENYLENE AND 
STILBENOID DENDRIMERS IN SOLUTION 

A. Form factors 

It has been demonstrated that the building units of 
the dendritic scaffold of polyphenylene dendrimers of the 
fourth generation (see Fig.[T](a)) are rather well-localized 
and no back folding of the terminal groups occurs'''^. A 
general feature of the polyphenylene building block is its 
finite angle formed by two subsequent phenyl group a^'^*^ . 
This is due to the strong repulsion between their o- 
hydrogen atoms that prevent a co-planar conformation. 
As a consequence, the polyphenylene dendrimers exhibit 
a three-dimensional structure which has been found by 
SANSl^. 

Such steric hindrance does not exists in dendritic scaf- 



P{q) = A dx 



— x'^) sin {qLx/2) 
[qRVT^f iqLx/2f 



(7) 



where Ji(x) denotes the cylindrical Bessel function of 
first-orderJ^. Equation ([7]) follows from Eq. ([5]) by re- 
placing the double sums by double integrals correspond- 
ing to a continuous distribution of interaction sites (n — > 
oo) and by using LUrnoiq, p) = sm{qlmo) / (qlmo) , where 
Uno — > is the bond length between sites m and o on the 
same cylindrical particle. 

The polyphenylene dendrimers were modelled by a 
central unit consisting of two small spheres (diameter: 
0.4 nm). This describes the central biphenyl unit (see 
Fig. [3]). The four dendrons were mimicked by eight 
spheres (diameter: 0.8 nm) with a center to center dis- 
tance of 0.87 nm. This distance is the approximate 
length of the two phenyl groups that connect two sub- 
sequent shells. Each sphere comprises five benzene rings 
that constitute a branching group. The endgroups are 
made of an equal number of phenyl rings. The form fac- 
tor of the polyphenylene dendrimers has been modelled 
by taking an average over approximately 500 conformers 
generated by randomly choosing the torsion angles for 
each dendron (see Ref4 for further details). The good 
agreement between the experimental and theoretical re- 
sults (Fig. (2) demonstrates that a full understanding of 
the spatial structure of both polyphenylene dendrimers 
of the fourth generation and stilbenoid dendrimers of the 
third generation has been achieved^ii^. The ability to 
construct coarse-grained models in such a way as to re- 
produce experimentally determined form factors is im- 
portant as it allows one to quickly study within PRISM 
theory various material systems based on input of a small 
amount of intramolecular information. 




FIG. 3: Modeling of the the form factor P{q) of the 
polyphenylene dendrimer of the fourth generation. The den- 
dritic scaffold is set up of a central unit consisting of two small 
spheres (distance: 0.4 nm) describing the central biphenyl 
unit. The four dendrons were modeled by 8 spheres (distance: 
0.87 nm), respectively, having a diameter of 0.8 nm. Both 
parameters were used as a fitting parameter. For the sake 
of clarity the size of the spheres has been depicted slightly 
smaller. The angles between the different units is indicated 
in the graph. The torsional angle was chosen at random for 
each dendron. 



B. Structure factors 



In Figs, m (a) and (b) the experimental structure fac- 
tors S{q, (f) are compared to the resuhs of the integral 
equation theory for the PRISM. We have used the cal- 
culated form factors P{q) (see the solid lines in Fig. ^ 
as input into the generalized Ornstcin-Zcrnike equation, 
i.e., P{qA) = P{q) in Eqs. §^ and (O. The particle 
number density is given by p — (j)/Vp, where is the 
volume fraction and Vp is the volume of an individual 
particle. The generalized Ornstein-Zernike equation is 
solved numerically together with the Percus-Yevick clo- 
sure. From Fig. |3]it is apparent that the PRISM integral 
equation theory is rather accurate. The magnitude and 
the scattering vector range of the suppression of S{q, p), 
i.e., the deviations from the value 1 at small scattering 
vectors, are characteristic for the size and the shape of 
the dendrimers as well as the volume fraction. The exper- 
imentally observed small upturns of S{q, p) at low scat- 
tering vectors q for polyphenylene dendrimers in solution 
(Fig. m (a)) indicate the presence of a small amount of 
aggregates due to attractive interactions which are not 
taken into account in the theoretical calculations. On 
the basis of our experience with both PRISM and RISM 
we expect that the results of the integral equation theory 
for the RISM would lead to very similar results, provided 
the same form factors are used^^. 
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FIG. 4: Experimentally determined structure factors for 
polyphenylene dendrimers of the fourth generation^ in (a) 
and stilbenoid dendrimers of the third generatioiji^ in (b) 
together with the results of the theoretical predictions of 
the PRISM integral equation theory (Eqs. Q and ((U). 
The volume fraction of the dendrimers increases from top 
to bottom: cf) = 0.016,0.032,0.046,0.065,0.08 in (a); (f) = 
0.009, 0.019, 0.039, 0.061, 0.1 in (b). For reasons of clarity, the 
upper data sets and lines have been shifted up by 0.1, 0.2, 
0.3, 0.4, respectively. 

IV. THERMODYNAMIC PROPERTIES 

The structure factor provides a direct link with ther- 
modynamics via the compressibility equatior>i^ 

lim S{q, p) = pknTKrip) , (8) 

where kt{p) is the isothermal compressibility. The os- 
motic pressure P{p) (equation of state) then follows from 

p 

^ = Jdp'S-\q^O,p'). (9) 


Various attempts have been made to develop accurate 
theories for the equation of state of fluids consisting of 
non-spherical particles: 

(a) Scaled particle theory^^iMi^^, which is very suc- 
cessful for hard sphere fluids, has been extended to 
prolate and oblate ellipsoids of revolution"^^ , however 
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with moderate success when gauged against Monte Carlo 
simulations'^^. Recently, it has been shown'^® that 
the results of scaled particle theor y^^'^^i"^° for plate- 
like particles or the closely related model of hard cut 
spheres are in disagreement with computer simulation 

(b) Onsager theory^i^S, based on the second viral co- 
efficient alone, can be " rescaled"^2i^i^. Also this semi- 
empirical procedure leads to reasonably good results for 
rodlike particles, it is much less satisfactory for platelike 
particlesS^. 

(c) Many theoretical studies on hard sphere fluids and 
depletion agents use the so-called free volume theory^, in 
which the free volume accessible to a single particle plays 
a major role. Recently, this free volume theory has been 
studied within a fundamental measure theory'^^. How- 
ever, it has been demonstrated that the resulting third 
virial coefficients of the equation of state for both hard 
cylinders and hard cut spheres differ from computer sim- 
ulation results (see tables I and II in Refi^). Theoret- 
ical approaches based on fundamental measure theory 
do not yield correct third virial coefficients and equation 
of states due to the occurrence of so-called " lost cases" , 
i.e., the fact configurations of three particles with pair- 
wise overlap but no triple overlap do not contribute to 
thermodynamic properties (see, e.g., Refsj ^^'^^ ). 

These earlier theoretical and computer simulation 
studies demonstrate that the understanding of thermo- 
dynamic properties of non-spherical particles needs to be 
improved. Here we model the measured inverse structure 
factor S~^{q = 0,p) extrapolated to vanishing scattering 
vectors of stilbenoid dendrimers of the third generation 
in terms of the so-called y3-expansion^2iS 

(10) 

« 1 + 2B2(f> + 3B3(t>'^ + 0{<P^) . (11) 

The 2/3 theory reproduces the exact second and third 
virial coefficients, B2 and B3, respectively. However, its 
practical applicability is limited due to the difficult nu- 
merical evaluation of the third virial coefficient B3 in the 
case of non-spherical particles^, while the second virial 
coefficient B2 for an isotropic hard convex body fluid is 
known exactly (see Ref.^° and references therein): 



B2 = ^{ V, 



ApRp 



(12) 



area and the mean radius, respectively, where the lo- 
cal mean curvature is denoted as Hp. For a circular 
platelet of radius R and thickness L one has Vp — ttR^L, 
Ap = 2'kR{R + L), and Rp = ttR/A + L/A. For platelike 
stilbenoid dendrimers of the third generation the second 
virial coefficient B2 — 5.54 as calculated from Eq. (fT^ 
with R = 2.4 nm and L — nm agrees with the exper- 
imentally determined second virial coefficient. 
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FIG. 5; Inverse structure factor S~^{q = 0, /o) extrapolated 
to vanishing scattering vectors of stilbenoid dendrimers of the 
third generationii (with the same symbol code as in Fig. |4] 
(b)). The dashed line follows from the scaled particle theory 
according to Eqs. ([10}, g!]), and ^ with Bf'"^'^ = 16.952 
while the dotted line represents the results of the fundamental 
measure theory as obtained from Eqs. (|10[) . (I12|l . and (|14|) 



with B 



(FMT) 



16.432. The solid line shows the results as 



obtained from Eqs. ^ and (HU with B3 = 34.129. 



In the framework of both scaled particle theory and 
Rosenfeld's fundamental measure theory S~^{q = 0,p) is 
also given by Eqs. (fTO| - (fT2|) but the third virial coeffi- 
cient is given by 



B. 



{SPT) 



1 



v. 



2RpApVp 



within scaled particle theory and 



B. 



{FMT) 



1 



-t- 2RpApVp -t- — 



(13) 



(14) 



within fundamental measure theory. In Fig. [5] 
the experimentally determined inverse structure factor 
S^^{q = 0,p) extrapolated to vanishing scattering vec- 
tors of stilbenoid dendrimers of the third generation is 
compared with the results of scaled particle theory ac- 
cording to Eqs. (fTUl) . (112), and ([T^ and fundamental 
measure theory according to Eqs. p0|) . ((T2|), and (fT4)) . 
With increasing volume fraction the theoretical results 
of both scaled particle theory (dashed line) and funda- 
mental measure theory (dotted line) deviate from the ex- 
perimental data (symbols). These deviations are mainly 
due to the fact that the predicted third virial coefficients 



Here Ap and Rp — (l/47r) J dAp Hp are the surface B 



(SPT) 



16.952 and B. 



(FMT) 



16.432 are too small. 



The results for S^^{q = 0, p) as obtained from Eqs. (fT0|) 
and (fT2|l with B3 — 34.129 as input are in agreement 
with the experimental data (see the solid line in Fig. [5]). 
We emphasize that the PRISM theory discussed above 
leads to a similar agreement with the experimental data 
as is apparent from Fig. [5] (b). 

Our first comparison of experimentally determined 
S~^{q = 0, p) of platelike particles with the predictions of 



6 



the well-known scaled particle and fundamental measure 
theory confirms earlier caveats concerning the applica- 
bility of these theories to freely rotating non-spherical 
particles. 

V. CONCLUSION 

We have presented a systematic application of the 
PRISM integral equation theory, scaled particle theory, 
and fundamental measure theory to rigid dendrimers in 
solution. The main findings and conclusions may be sum- 
marized as follows. 

PRISM theory is quite successful in reproducing ex- 
perimental structure factors S{q, p) of both polypheny- 
lene dendrimers of the fourth generation [Fig. 2] (a)] and 
stilbenoid dendrimers of the third generation [Fig.[4](b)], 
provided the correct form factor P{q) [Fig. [5] is used as 
input into the generalized Ornstein-Zernike equation of 



the PRISM [Eq. ([8])]. These investigations encourage to 
pursue a study of charged dendrimers^^'^^'^'^ within the 
framework of PRISM by using interaction sites which 
carry charges. 

The inverse structure factor S^^{q 0, p) extrapo- 
lated to vanishing scattering vectors as predicted by both 
scaled particle theory and fundamental measure theory 
is systematically smaller than the experimental data for 
platelike stilbenoid dendrimers of the third generation 
and the prediction of the PRISM integral equation the- 
ory [Fig.[5] . The substantial differences observed between 
the experimental data and the results of scaled particle 
theory and fundamental measure theory are mainly due 
the fact the both theories do not yield the correct third 
virial coefficient in the case of platelike particles. Hence 
there is a clear need to improve both scaled particle the- 
ory and fundamental measure theory for freely rotating 
non-spherical particles. 
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